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Abstract 

We formulate a linear difference equation which yields averaged semi-inclusive 
decay rates for arbitrary, not necessarily large, values of the masses. We 
show that the rates for decays M — > m + M 1 of typical heavy open strings 
are independent of the masses M and m, and compute the "mass deffect" 
M—m—M'. For closed strings we find decay rates proportional to Mm^ D " 2 ^ 
where m,R is the reduced mass of the decay products. Our method yields exact 
interaction rates valid for all mass ranges and may provide a fully microscopic 
basis, not limited to the long string approximation, for the interactions in the 
Boltzmann equation approach to hot string gases. 
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I. INTRODUCTION 



The decay properties of very massive strings were actively studied in the past 
and several general features, such as the existence of a constant decay rate per unit length 
and the dependence of the total decay on the mass of the string emerged. Many of those 
works extracted the decay rates from the imaginary part of the self-energy diagram and, 
for technical reasons, the computations were restricted to decays of very particular states, 
such as strings on the leading Regge trajectory |T]-^[ or strings stretched around a compact 
spatial dimension ||. Moreover, given that the one-loop diagram involves an integral over 
the momenta of the decay products, no information was obtained on the energy distribution 
of the emitted strings. 

However, we can easily imagine situations where we don't know or are not interested in 
the full details of the initial state. For instance, a high temperature string gas will contain 
strings of many different masses, and we may want to find the corresponding distribution 
function |^-[7j . One possible approach is to try to formulate and solve a Boltzmann equation 
for the system PH, and for that we need to know the decay (and recombination) rates of 
typical strings of mass M. The averaged total interaction rate of two highly excited closed 



strings was computed in |10[ and, for massive long strings in different situations, one may 
usually guess the form of the rates from semiclassical considerations But it would be 

interesting to have a fully microscopic formulation valid for all mass ranges. 

Similarly, a highly excited string is believed to form a black hole at strong coupling. But 
we can try to mimic this at weak coupling by pretending that we only know the mass of 
the state, and see if the properties obtained by taking averages over all the states of given 
mass M resemble those of an actual black hole. In particular, we would like to compare the 
energy and mass spectrum of the emitted strings. 
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In this spirit, the photon emission rate by (bosonic) massive open strings^ was considered 

in 



|T2|j . This can be formulated as an initial string state going into a final state |$/) 
through the emission of a photon of momentum k 11 and, at weak coupling, the amplitude 
for this process is proportional to the matrix element (§f\V(k)\&i) of the photon vertex 
operator. Q The resulting decay rate, which can be computed for any initial and final states, 



has in general a complicated polynomial dependence on k^. However, the authors of fi2 
showed that, when these rates are averaged over the many initial states with mass M 
and fixed momentum p M and summed over all final states [$/), the resulting spectrum is 
thermal. More precisely, they considered the averaged semi-inclusive emission rate 

1 



i v 



dT N (k°) « — EE \(^ f \V(k)\^)\ 2 k D - 3 dk , (1.1) 



where the sums run over all states satisfying the mass conditions a'M = N — 1, 
a'M' 2 — N' — 1, Q{N) is the degeneracy of the N th mass level, and D is the space-time 
dimension. Evaluating ( |1 . 1| ) for large N they obtained the following result for the photon 
emission rate 

_j£_ 

dT N (k°) = const. 6 k D ~ 2 dk , (1.2) 

1 - 

which is the spectrum of a black body at the Hagedorn temperature 



T H =—^\— — ■ (1-3) 



This identifies the Hagedorn temperature [15] as the radiation temperature of an averaged 



string, and can be seen as an example of the emergence of classical or statistical concepts 



from microscopic computations through the average (decoherence) process [12.16]. It also 



1 They also studied massless emission by charged closed strings and obtained grey body factors of 
the type found in black holes [11]. We will consider closed strings in Section 4. 
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See, for instance [1^,14[. 
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shows that the properties of typical or averaged strings can be very different from those 
exhibited by specific states. 

It is somewhat surprising that all previous studies regarding the decays of fundamental 



strings missed this result, but the calculation in |12| differs in two important points from 
those carried out in the past [[]]-[§] . Namely, rather than fixing the initial state and summing 
over all final decay products, (|1 . 1|) fixes the state of one of the emitted strings, and sums over 
all initial and final states of masses M and M'. Also, unlike previous computations (|1 . 1| ) does 
not sum over the momenta of the decay products, thus yielding the energy spectrum of the 
emitted photon. If we can evaluate ( |1 . 1|) not only for photons, but for any vertex operator 
V(k), we will have complete information on the emisssion spectrum of typical strings. In 
particular we will be able to compare the rates at which strings of different masses and spins 
(representations) are emitted. 

Having such a detailed information on the decay modes of typical strings can be useful 
in several ways. For instance, we may try too see which of the features present in the 
decays of strings on the leading Regge trajectory [|TJ-g] apply also to typical strings. This is 
important, since states on the leading Regge trajectory are very atypical, maximum angular 
momentum states. Other possible applications include a fully microscopical derivation of 
the interactions appearing in the Boltzmann equation approach to string gases HH, and 
a comparison between the decay properties of strings and black holes. For masses below 
the correspondence point \ Tf\ , black holes transition to a random- walk-size string state or 



"string ball" |T8|,[T9|, which may even be accessible to production in the LHC or VLHC |T9 
If this turns out to be the case, computing its decay properties is obviously relevant. 

These are the main motivations behind this paper, where we develop a method which 
allows the computation of the double sum in ( |1 . 1| ) for any (not necessarily large) values of 
N and N', and for vertex operators corresponding to arbitrary physical statesQ As we will 



For simplicity, only the bosonic string is considered in this paper. 
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see, a direct evaluation of becomes very involved when V(k) corresponds to a massive 
string, but the problem can be reformulated in terms of a linear difference equation. This 
equation depends only on the mass of the emitted states, and different "initial conditions" 
correspond to different states with the same mass. In this way we can neatly separate the 
universal properties of the decays from those dependent on the details of the emitted state. 
In particular, we can show that as long as the mass m of the emitted string is less than 
2M/3 (M is the mass of the decaying string), the emission rate is essentially independent of 
its state. The averaging process leaves only a universal dependence on m. This is no longer 
true for m > 2M/3. 

The solutions to the linear difference equation can be used to study the decays of open 
and closed strings, both heavy and light. For example, we show that the rates for decays 
M — > m + M' of heavy open strings are independent of the masses M and m, and compute 
the mean value for the "mass deflect" M — m — M' . For closed strings we find decay rates 
proportional to Mm|| D )/ 2 ) w here uir = m(M — m)/M. Since our solutions are valid also 
for small values of M, we can see how the decay properties evolve as M changes. 

This paper is organized as follows. In Section 2 we present a brief review of photon 
emission as computed in |T2| and obtain a general formula which applies to the emission of 



arbitrary massive states. This formula involves contour integrals over a modified correlator 
on the cylinder, and a direct evaluation is difficult . However, in Section 3 we show that the 
integrals satisfy a universal recursion relation, which depends only on the mass of the emitted 
particle, and compute and analyze the solutions. These solutions are used in Section 4 
to study the processes of fragmentation and radiation by heavy open and closed typical 
strings, and a comparison is made with previous results for strings on the leading Regge 
trajectory. We also analyze decays of (relatively) light strings. Our conclusions and outlook 
are presented in Section 5. 

II. MASTER FORMULA FOR EMISSION OF ARBITRARY STATES 
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A. Review of photon emission 

The first step in the evaluation of the averaged inclusive rate in jT2] is the transformation 
of into a double complex integral. To this end, write the double sum as 

F(N,N') = £ E \(9,\V{k)\9 t )\ a = E E (^(k)\^ f )(^ f \V(k)\^) (2.1) 
and introduce projection operators over the mass levels 

P N=( f^ z N-N AT =E (2.2) 

With their help the sums can be converted to sums over all physical states, yielding 

F(N, N') = I —z~ N I ^-z'- N 'Tt [z*V\k, l)z'*V(k, 1)] , (2.3) 



w z Jc> z' 

where C and C are small contours around the origin. The vertex operator has been arbitrar- 
ily placed at z — 1, although SL(2, R) invariance of open string tree amplitudes guarantees 
that the result is independent of this choice. For photons 

V(k,z) = i^d T X^{z)e ik - x{z \ ^k" = 0, k 2 = , (2.4) 

where we have introduced the Euclidean proper time r = In z and X M is given by the mode 
expansion?] 

00 1 

X»(z) = x» - ijf In z + % E - (<z~ n ~ a-n^j = x" - ip" \nz + X£{z) 

n=l n 

\x»,lT\=iTr , K, «_J = nWT ■ (2-5) 



Notice that we have defined the oscillator part Xg(z) for latter convenience. The trace 
in (|2.3|) is evaluated by coherent states techniques and, making the change of variables 



v = z' and w = zz' , the result isf| 



4 Henceforth we set a' = 1/2. 

5 In order to obtain this result one should also introduce coherent states for the ghost variables 
and add the corresponding contribution to the number operator in the projectors (|2.2j ). 
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F{N, N') = U» f —w~ N f(wf- D I -v N ~ N ' ( Pfl p„ + Vlu ,n(v, w)) , (2.6) 

JC W W JCv v 

where 

and f(w) is related to the Dedeking //-function by 

f(w) = mi-w n )=w- 1 ^(^) ■ (2.8) 

n=l V Z7FZ / 

Note also that the power of f(w) in fl2.6| ) reproduces the partition function for the open 
string, and its power expansion generates the mass levels degeneracies 



Ttw n =U(l-w n ) 2 ~ D =J2S(N)w N . (2.9) 

n=l N=0 



The f-integral can be done trivially, and we are left with 



F{N , N < )=e(N - N <) ^ N ^r D . ( 2 ,o) 



For large N' this integral can be computed by a saddle point approximation, with the 
main contribution coming from w ~ 1, where f(w) may be approximated by 



_2 



f(w) ~ const. (1 - w)~ L/2 e 6 ^~ w ' . (2.11) 
This, together with the asymptotic formula for the degeneracy of mass level N 



D+l 



'D-2 



G(N) ~ N~^e aVW , a = 2n ] l^-^, (2.12) 



leads to the black body spectrum fll.2|) mentioned in the Introduction. 

The computation above considers only one of the two cyclically inequivalent contributions 
to the three-point amplitude. However, these two contributions are related by world-sheet 
parity, and differ only by a sign equal to the product of the world-sheet parities of the states. 
For the open bosonic string the parity of a state depends only on its mass level and is given 
by (— ) N . Thus for photon emission the relative sign between the two cyclically inequivalent 



contributions is given by (— ) iV + 7V '+ 1 and, for £7(1) strings, the total amplitude is zero for 
N + N' even. For other groups the decay amplitude for iV — > Ni + N2 should be multiplied 
by tr[A(AiA2 + (— ) < - Ar+Ari+7V2 - ) A2Ai)]. For the rest of the paper we will assume a U(l) group, 
and keep in mind that decay rates vanish for N + Ni + N 2 odd. 

Before considering the emission of arbitrary states, we would like to make a remark re- 
garding the meaning of (|1.1| ) for D 7^ 26. Since we are not considering a compactification, 
strictly speaking the theory makes sense only for D = 26. However, this is just a tree 
level computation, and it seems that as long as we restrict ourselves to the "transverse", 
positive norm states generated by D — 2 oscillators we should not find any inconsistency. 
This restriction is usually described as a "truncation" of the theory, and was widely used 
in previous computations of decay rates [I]]. At the critical dimension, the total number 
of propagating physical states f\ coincides with those generated by the 24 transverse oscil- 
lators of the light-cone gauge. However, for D < 26 there are some extra "longitudinal" 
propagating states which can not be obtained from D — 2 oscillators, and ( [2.9|) does not 
give the total degeneracy, only the number of "transverse" states^]. This would not be a 
problem if the sums in fl2.1|) counted only transverse states, as the power of f(w) in ( |2.6|) 
would suggest. However, a explicit computation in Section 4 shows that at least some of the 
longitudinal states are counted in F(N,N'). On the other hand, the recursion relation to 
be obtained in Section 3 applies only to the emission of transverse states. Thus, it seeems 
that one can have an entirely consistent picture only at the critical dimension, even though 
only tree amplitudes are beeing computed. Although we will continue to write most of our 
formulas for arbitrary D, one should bear in mind that their validity is clear only for the 
critical dimension. 



6 By physical state we mean one that satisfies the Virasoro conditions (L n — 6 n o)\$) = 0, Vn > 0. 
Propagating states are physical states with positive norm. 
7 The problem of counting physical states for D < 26 is expained in great detail in [20|. 
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B. Emission of arbitrary states 



As is well known, to any physical state 

l*> = Ew»«V--«XI ;*> ( 2 - 13 ) 

we can associate a vertex operator V(k, z) by making the following substitutions [14 

iff^y (o^^^w (2 . 14) 

(n — 1)! 

and normal ordering the result. The simplest such state is the tachyon 

V(k, z) =: e ik - x(z) : , A; 2 = -m 2 = 2 . (2.15) 

Since Q2.3|) is valid for any physical state, provided that we substitute the appropri- 
ate vertex operator, we can use the same standard coherent states techniques for tachyon 
emission. Instead of (E.6I) we obtain 



F(N,N') = f—w- N f(w) 2 - D l—v N - N ' ^(v,w)- 2 (2.16) 

J W J V 

where 

iftv, w )= 1 - v I] " n ~ w ^ 2 ' ' ■ 2.17 

„=i (1 - w™) 2 

Obviously, the form of the function ip(v, w) renders a direct evaluation of the w-integral 
rather non-trivial. This point will be addressed in the next Section. Here we will consider 
the problem of writing the equivalent of ( |2.16| ) for an arbitrary physical state in a systematic 
way. One key observation is that tp(v, w) is related to the scalar correlator on the cylinder 

(X^l)X u {v)) = -tT^(v,w) (2.18) 

by 

f-\n 2 v\ 

ip(v,w) = y/vexp — ) ip(v,w) . (2.19) 

\ 2 In w J 
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Similarly, we recognize that the integrand in (|2.16|) would correspond to a two-point ampli- 



tude on the cylinder if we made the substitution 



j^r D ^f(wr D (-^-) D/2 . (2.20) 



Into. 

These are not accidents. In the operator formalism] 5 ] the integrand for the two-point 
one-loop amplitude is given by 

d D pTr [z L °V\k, l)z °V(k, 1)] = f{wf- D [ ) (VUk, 1) V(k, v)) , (2.21) 

\mw J 



where v — z' and w = zz' . The extra factors in ( 2.19 ) and (|2.20|) are the result of the 



momentum integral over the zero mode parts of the Virasoro generator L , which do not 
appear in the definition (|2.3|) of F(N,N'). This implies that, for any physical state, the 
evaluation of the trace in F(N, N') will give 

Tr [z*\rt(k,l)z!*V(k,l)]=Tr [V'\k, l)V'(k, v)w*] 

= f(w) 2 - D (V'\k, 1) V'(k,v)) , (2.22) 

where the prime over the vertex operators means that the following substitution should be 
made 

d?X»(v) = (xP - ip» In v)5 n , - ijf S n ,i + d n T X»{v) 

-> v*d?X>*(l)v-" = 7? 5 nfi - ip" 5 nA + d?X£(v) . (2.23) 

The reason for this substitution is that, although Lq is the generator of proper time trans- 
lations 

z Lo V(k,l)z- Lo =V(k,z) , L = ^p 2 + N, (2.24) 

the operator N does not affect the zero modes in d^X^{l). Note that the prime in V'(k, 1) 
is irrelevant. 



8 See for instance chapter VIII of [ 13 1 
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The correlator in (|2.22|) can be computed by using 



(d?Xg(l)d?XZ(v)) = -(-) n r]^(vd v ) n+m \niP(v,w) 



and the master formula for emission of general states is simply given by 



(2.25) 



r dw 



dv 



F(N,N') = & —w~ N f(w) 2 ~ D & —v N - N ' (V'Hk,l) V(k,v)) 

J W J V 



(2.26) 



which should be evaluated according to the prescriptions in ( |2.23| ) and ( |2.25| ). It is easy to 
check that the tachyon fl2.16| ) and photon fl2.6| ) sums follow immediately from this formula. 
One should simply notice that 



/ e -ik.X (l) e ik.X (v)\ _ e ~k 2 \ntj)(v,w) 



ip(v, w) 



(2.27) 



and 



Q(v, w) = — (vd v ) 2 \nip(v, w) 



(2.28) 



We shall illustrate the use of the prescriptions in fl2.23|) and (|2.25|) with two further 
examples. The vertex operator for the symmetric rank-two tensor at mass level N = 2 is 
given by 



V(k, z) = ^ v : d T X»(z)d T X u (z)e ik - x{z) : , k 2 = -m 2 = -2 

1 

2 ' 



(2.29) 



where the last condition on the polarization ^ u guarantees that the state is normalized. We 
obtain 



(V'\k,l) V'(k,v)) = ^P(v,w) 2 

(W) 2 + 4(£ M pC WW*, w) + 2«^>, w) 



x 



(2.30) 



Similarly, for the rank-two antisymmetric tensor at mass level N = 3 



V(k,z)=^ u :d T X»(z)d 2 X»(z)e lk - x ^: , k 2 = -m 2 



(2.31) 
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the result is 



(V'\k,l) V'(k,v)) = i){v,wf 

^pC W«(*>, w) + ^rmv, w)Q(v, w) - (l{v, w) 2 ) 



x 



(2.32) 



where the dots stand for derivatives with respect to proper time (d T = vd v ). Obviously a 
direct evaluation of the v-integral for these correlators would be very difficult. 

We close this section by pointing out that more complicated vertex operators may require 
evaluations of self-contractions between operators at coincident points. These are easily 
computed by using 



where 



-) n Y u (vd v ) n+m \n^ R (v,w) 



tp(v, w) 



v=l 



(2.33) 



[l-v) 



(2.34) 



Note that the denominator cancels the zero at v — 1 and all the derivatives appearing 



in (|2.33|) are non-singular. This simply reflects the fact that vertex operators are normal 
ordered. 



III. RECURSION RELATION 



A. The derivation 



In this section we address the problem of computing the ^-integral in the master for- 
mula (|2.26| ). We will show that, if we rewrite the master formula in terms of new functions 



IJw) 



l n (w)=I -v n (V'\k,l)V'(k,v)) 
Jc v V 

dw 



F(N, N') = I —w~ N f{w) 2 - D I N „ N ,{w) , (3.1) 
Jc w w 
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then the properties of the correlator imply a of recursion relation for I n (w). In order to 
simplify the derivation, we will write the vertex operators in a very specific gauge. Namely, 
the coefficients c Ail ... Atm in ( |2.13| ) will be orthogonal to the momentum k^ 



Note that the vertex operators (|2.4j) , (|2.29| ) and (|2.31|) are written in this gauge. Although 



for photons (|3.2j ) is a consequence of the physical state conditions, it becomes a gauge choice 



for states with N > 2. This gauge was analyzed in detail in pl| , where it was shown that for 
each propagating (i.e. positive norm) physical state there is a unique gauge representative 
satisfying (|3.2| ) for D = 26. For D < 26 there are positive norm, physical "longitudinal" 
states which can not be written in this gauge, and our formalism is not useful in computing 
their emission rates. 

For vertex operators satisfying ( |3.2| ) the only contraction involving the exponents is given 
by eq. ( |2.27| ). All other contractions will involve at least two proper time derivatives, and 



the correlator in X n {w) will take the following form 

(V'\k, 1) V'(k, v)) = ij m2 V(Q, Q, Q, . . .) , (3.3) 



where V is a polynomial (cf. (gj), (gH ), (|2~30|) and (ggj) ). 



On the other hand, the contours in (|3.1|) must satisfy |iw| < \v\ < 1, which follows from 
\z\, | z' | < 1 and v — z', w — zz' . The functions T n (w) can be viewed as the coefficients in 
the Laurent expansion of the correlator 



(V'\k,l)V'(k,v))= J2 U^)v n (3.4) 

and this expansion is valid only inside an annulus between the nearest singularities, which 
are the poles at v — 1 and v = w.f\ 



9 A general account of this and other properties of the correlator, as well as its relation to Jacobi 
theta functions can be found in Appendix A. 
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FIGURES 




FIG. 1. Contours used in the derivation of the recursion relation (I 



Let us consider a second v -contour C' v denned by 1 < \v'\ < \w\ 1 , which is related to C v 
by the transformation v' = w~ 1 v. Using the properties 



A (J Ss 

tl)(w~ v,w) = ip(v,w) 



w 



d?tt(w- l v,w) = d?n{v,w) 



(3.5) 



which follow directly from the definitions of ip(v,w) and Q(v,w) and taking into account 
the general form of the correlator ( p.3|) , we may relate the two contour integrals 



dv 



C> V 



v n {V'\k,l) V'(k,v)) =w 



dv 



1 " -v n+m2 (V'\k,l) V'(k,v)) 

ICv v 



w 



-n—m 



■T 

■ J -n+rn? 



» • 



(3.6) 



The difference between the two contours can be deformed into a small contour C around 
v — 1 (see Fig. 1), and we have 



dv 



—v n (V'\k,l) V'(k,v)) =w~ 
C' v Jc v l v 



l n+m 2(w) -l n {w) 

—v n (V'\k,l) V'(k,v)) . 

C V 



(3.7) 
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This last integral depends only on the structure of the singularity at v = 1 , and is computed 
in Appendix A. Assuming that the emitted state is normalized, the result is 

/ — v n (y'i(k, 1) V\k, v))=n+ T ^ r (3.8) 

and, after multiplying ( p.7| ) by w n+m " 2 , the following recursion relation is obtained 

2 

9 77? 

l n+m2 (w) = w n+m [X n (w) +n + —]. (3.9) 

This relation is universal, in the sense that it depends only on the mass of the emitted 
state. By a similar argument involving two contours related by the transformation v' = w/v 
the following "reflection symmetry" is easily derived 

X n (w) = w n X_ n {w) . (3.10) 



B. The solutions 

For photons, the recursion relation reduces to X n {w) = w n [X n (w) + n] with the obvious 
solution 



X n (w) = , (3.11) 

1 — w n 



which agrees with (|2.10 ) and gives rise to the Plank distribution. For m ^ there are 



no closed-form solutions, but they can be expressed as power power series in w. Since 
"4>(v, w) and all proper time derivatives of Q(v, w) have finite limits for w — > 0, only positive 
powers of w are involved. The absence of negative powers of w combined with the reflection 
symmetry ( |3.10| ) implies 



X n (w) = a n w n + 0(w n+l ) (3.12) 

One can easily check that, for tachyons (m 2 = —2), this condition implies that the unique 
solution to (|3.9p is given by the series 
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Z n {w) = J2( n - 1 + 2p)w {n ~ 1)p+p2 . (3.13) 
P =i 



For positive m 2 = 2, 4, 6, ... we can use (|3.9| ) to express X n {w) in terms of X n „ m 2 (w), which 
in turn can be related to X n _ 2m 2(w) and so on, until one reaches X v [w), where < v < m 2 . 

The result of this iterative process is the following series 

a 1 

X n ( W ) = E - m2 (P - -))w np -^ m2{p2 ~ p) +1 U {W) w vA+\m\K>+A) ^ (3 W ) 

p=l 2 

where v = n mod m 2 and n = m 2 A + v. At first sight the solution ( |3.14|) depends on the 
m 2 undetermined functions T , Xi,---X m 2_ 1 . However, combining the reflection symme- 
try ( |3.10| ) with the recursion relation yields 



2 

I u {w) = w 2v ~ m l m 2_ u (w) + (u-—)w\ (3.15) 

which implies that only X , X 1; • • •X m 2 / / 2 are independent functions. 

Note that the structure of this solution is consistent with the fact that ( p.9| ) is an in- 
homogeneous linear difference equation. The last piece in ( |3.14j ) is the general solution to 
the homogeneous equation X n+m 2(w) = w n+m,2 X n (w), which exists only for m 2 > 0; the sum 
from p = 1 to A is a particular solution to the complete equation. 

As mentioned above, the recursion relation is universal, in the sense that only the 
mass of the emitted particle enters it. However, the solutions depend on the functions 
X , Xi, • ■ -X m 2/ 2 , which should be determined by direct expansion of the correlator. Since 
this expansion depends on the particular vertex operator considered, we may say that the 
solution to the homogeneous equation contains the non-universal, state dependent part of the 
emission rate, while the particular solution reflects the universal behavior. The non-universal 
terms exist only for m 2 > 0. This is physically sensible, since there is only one tachyonic 
state and all D — 2 transverse physical photon polarizations are equivalent by SO(D — 2) 
invariance. 

In order to compare the relative strength of universal and non-universal contributions to 
the decay rate, we should note that if X n is given by its power expansion X n {w) = J2i>o a i w \ 
then by (p|) and Q F(N, N') will be given by 
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F(N,N')= y £a l g(N-l) . (3.16) 

l>0 

Since the degeneracy Q(N) is a rapidly growing function of N, the higher powers of w will be 
strongly suppressed. On the other hand, the highest power in the universal series in (|3.14|) 
is the term with p = A, which gives / = vA + ^m 2 (A 2 + A). This is precisely the power 
of w which multiplies X v {vS) in the non- universal term and since T v (w) = 0(w u ), we see 
that the non-universal term will tend to be suppressed in relation to the leading universal 
contributions. This argument fails only for A = which, as we will see in the next Section, 
corresponds to the emission of states carrying a large fraction of the total mass M. 



IV. APPLICATIONS 

In this Section we will describe some applications of the solutions to the recursion rela- 
tions. Substitution of the solution (|3.14j) into ( p.l|) gives the sum over initial and final states 



as a linear combination of level degeneracies 

A 1 1 

F(N,N') = Y J {n-m\p- )\q\n - np + -m\p 2 - p)] +F NU (N,N') , (4.1) 

p=i ^ ^ 

where n = N — N' and F^u is the contribution from the non-universal term 

F NU (N,N') = £ —w- N f{wf- D XJw)w vA +^ A ^ . (4.2) 
Jc w w 



Note that ( |4.1|) is also valid for photons and tachyons if we set A = oo and Fnu — 0. The 
mass level degeneracies can be obtained by direct expansion of the partition function ( |2.9| ) 
or from an approximate formula, given in Appendix B, which is more accurate than the 
usual asymptotic expression fl2.12p . 



Taking into account the normalization of the initial state and the two-body phase space, 
the averaged semi-inclusive decay rate is given by 

r(M - m + M') = Ao^^j^ k D ~* , (4.3) 

where 
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and g a is the open string coupling constant. Remember that the label "averaged semi- 
inclusive" refers to the fact that ( |4.3| ) describes the production of strings of mass m in a 
very specific state, by typical strings of mass M, while summing over all final states of 
mass M' . 



A. Fragmentation of heavy open strings 

Here we will assume that the mass of the initial string is so large that we can use the 
following asymptotic approximation for the mass level degeneracies 



G{M) ~ (2nT H )- ^M-^eWT* , T H = - J . (4.5) 

it \ D — 2 



This is ( j2.12 ) written in terms of the mass, and we have specified the value of the multi- 



plicative constant in order to get quantitative results. Using ( fOD and 

n = N-N' = Mk°- ^m 2 > m(M - ~m) , (4.6) 
which follows from p = p' + k and the mass shell conditions, we can obtain the following 



estimate for the ratio between the first two terms in flO 



g(N -2n + m 



2\ 



exp 



m /M-3m/2\ 
\ M-m J 



(4.7) 



G(N-n) 

Obviously, as long as < m < 2M/3 the contribution of all the terms but the first can be 
neglected, and we may take 

F(N, N') &(N-N' - im 2 ) Q(N') . (4.8) 

If we further asssume that m <ti M we can use the approximation n = N — N' ~ Mk° 
together with Q{N')/Q{N) w e~ k ° ITH . Then (fO|) implies the energy distribution 

dT(k°) w A D g 2 e- k ° /TH k D - 2 dk , (4.9) 



where we have used dn = Mdk° and k°dk° = kdk. This confirms that massive strings 
with m <C M are radiated according to a Maxwell-Boltzmann distribution which is totally 
independent of the mass M of the decaying string. Had we included the rest of the terms in 
the series in F(N, N') we would have obtained Plank distribution. But even for the lowest 
lying massive states k° > m is large compared with Tjj, and the difference between the two 
distributions can be neglected. 

In what follows we are going to lift the restriction m <C M and consider instead the 
splitting of a string of mass M into arbitrary states of masses m and M'. Note that ( |4.3|) 
gives the average rate of decay of a string of mass M into a specific state of mass m. However, 
our approximation (|4.8|) for F(N, N') implies that all the states of a given mass level m are 
emitted at the same rate. Thus, if we multiply the rate by Q (m) the result will be the total 
rate for decays into arbitrary states with masses m and M' 

V M (M - m + M') « A D -§- 2 (N - N' - ^-)^^ . (4.10) 

Note that, unlike (|4.8|) , this formula is valid even for m > 2M/3. The reason is that ( 4.10[) 
is symmetric under interchange of the emitted strings, and at least one of the two emitted 
masses is going to be less than 2M/3. The only conditions for the validity of ( f4.10|) are that 
M sould be large and neither m nor M' should be zero. 
As long as m > 1 and M' ^> 1, we can use (|4.5|) to write 



gim)g(M') . m , D-i mM'\ 2 t/r , 
K \ ' ~ {2nT H ) — - e - t/TH , 4.11 

g(M) \ m J ' v ; 

where t = M — m — M 1 is the total kinetic energy or "mass deffect" for the decay. Since 
the exponential dependence on t implies t ~ 0(Th) M, we can also make the following 
approximations 

N-N' w mM' , t « - f ) k 2 (4.12) 

and find for the total decay rate 

T tot (M ^m + M')^ A D g 2 { ^r~f t ^ e -t/T H (4 13) 

2mM' 
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Thus we obtain a Maxwell-Boltzmann distribution for the total kinetic energy. Note that 
if we only measure the energy of one of the decay products (say the one with mass m) we 
will find a temperature 

r„, = r H (i-B). (4 . 14 ) 

Of course, this is simply a "recoil effect", and in the limit m « M we recover our previous 
result (|4.9|) . Another consequence of the distribution ( 4.13|) is that the average kinetic energy 
released in a decay is given by 

< t >= ~(D - l)T H . (4.15) 

In other words, as long as neither of the final masses is too small, the energy released in the 
decay M — > m + M' is independent of the masses and satisfies the equipartition principle in 
D — 1 spatial dimensions. 

Another quantity that we may compute from ( |4.13| ) is the total rate of production of 



strings of mass m. We simply have to sum over all values of N' using dN' = —M'dt, and 
we get0 

; (7rT^)-^ r°° ,_d-5 _ t/2 , A' D g 2 Q 



m 



E r to ,(M -> m + M') ^ A D gl ^ lH J * f°° dtt^e^ = ™ , (4.16) 
w 2m Jo m 



where A' D = n 2 ~ D 2~ ( ' D+1 \ In terms of the level density p(m) — m, which follows from the 
mass shell condition, there are p(m) dm levels with masses between m and m + dm. Thus 
the number of particles produced within that interval is given by the density 

^^ = T m p(m)^A' D gl, (4-17) 

which is independent of the masses. This means that a heavy open string radiates strings 
of all masses with equal probabilities and implies a total decay rate proportional to the 



10 Since a bosonic string can decay into heavier states throught tachyon emission, the total rate 
diverges. Following |l||3| we avoid the problem by imposing N' < N by hand, i.e. by restricting 
ourselves to non-tachyonic decays. 
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mass of the string. Thus, the decay of a heavy open string can be described as a process of 
fragmentation, rather than one of radiation. 

From a semiclassical point of view, this can be explained in terms of a constant splitting 
probability per unit of length |3|,|]]. Since the released energy ( |4.15| ) is small compared to 
the masses, M « m + M' and all masses smaller than M are produced with the same 
probability. This should fail when one of the masses is comparable to the released energy. 
From the point of view of our computation, this failure can be attributed to the breakdown 
of the asymptotic approximation for the level multiplicity. 

The total decay rate of a string of mass M is obtained by integrating ( [4.1 7|) from to 



M/2, to avoid counting each decay mode twice. The result for D = 26 is g 2 ir~ 2A 2~ 28 M 
and can be compared to other computations in the literature. The decay rate for open 
bosonic strings on the leading trajectory is given by eq. (12) of [Q] and exactly agrees with 
ours. They obtain their result from the one-loop planar diagram, which is equivalent to 
considering only one cyclic ordering for the three-point tree amplitude. Including the non- 
planar contribution as in |T2] amounts to adding the other cyclic ordering, but the result 



is easy to predict. As mentioned in Section 2, amplitudes for U(l) strings vanish or get a 
factor of 2 depending on the parity of N + Ni + N 2 . The net result is that the total decay 
rate should be multiplied by 2, yielding g 2 n~ 2i 2~ 27 M . This can also be compared with the 
decay rate for a long straight open string^ , which according to eq. (14) of |4j] is given by 
g 2 7r~ 24 2~ 26 M . This is twice our result, but the discrepancy may be due to some overlooked 
subtlety. In any case, since we are computing the averaged decay rate for strings of mass 
M, the result need not coincide with those obtained for highly atypical states such as those 
on the leading Regge trajectory. 

We can also use ( |4. 1Q ) to compute the production rate T m by strings of mass M numeri- 



cally and compare the results with the asymptotic prediction ( |4.17| ). Fig. 2 shows our results 



ii 



Our oriented £7(1) string should be equivalent to the unoriented SO(2) string considered in 
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for N = 20 and N = 100 in 26 dimensions, with dV j m > normalized so that 1 corresponds to 

am c 

the asymptotic value ( |4.17|) . Note that even for N as low as 20 the rate manages to attain 
the predicted value for m ~ 2.3 and that for N = 100 there is already a well defined plateau, 
with a decrease in the rate near the ends of the mass range. 



0.8 
. 6 
0.4- 
0.2 



1 
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0.2 



10 12 14 



FIG. 2. for D = 26 > N = 20 ( M ~ 6 - 16 ' left ) and N = 100 ( M ~ 14 -!> ri e ht )- 



For D 26 the results do not match ( |4.17| ) so well. For D < 26 the production of light 
strings is enhanced, while the opposite is true for D > 26. This is illustrated in Fig. 3, 
which suggests that strings decay according to a the constant splitting probability only in 
26 dimensions ||. However, one should not forget that the consistency of the formalism is 
suspect for D ^ 26. Since Q(N) does not count the total number of propagating states for 
D < 26, and for D > 26 there are also negative norm physical states, the validity of ([4.1 0|) 
for D ^ 26 is unclear. 
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FIG. 3. for N = 500 and D = 16 (left), D = 26 (center) and D = 30 (right). 
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B. Radiation by heavy closed strings 



Our previous computations extend easily to the case of heavy closed strings but, as we 
shall see, the results are remarkably different. Tree amplitudes for closed strings factorize 



into a product of left and right open string amplitudes , and instead of (f4.3|) we have 



T C (M - m + M>) = A D jLE<W¥ , (4.18) 

where g c is the closed string coupling constant, which is proportional to g 2 . For closed 
strings, the mass shell condition becomes a'M 2 = 4(N — 1), and this introduces numerical 
factors into some of the previous formulas. For instance, eq. (|4.6|) becomes^ 

n = N-N' = ^{Mk° - ^m 2 ) . (4.19) 

The computations proceed along the same lines and we will merely quote the main 
results. For m <C M the energy distribution of emitted strings is given by 

dT c (k°) « -^A D g 2 c M k°e- k ° /TH k D ~ 2 dk . (4.20) 

Besides the trivial numerical factor, there are two important differences with the analogous 
formula ( [4.9|) for open strings. These are the presence of the mass M and the factor A; . The 
latter can be seen as a grey body correction to black body spectrum, similar to the ones 
found in other contexts [fLl]| . The factor of M indicates that, unlike open strings, closed 
strings radiate proportionally to their mass. 

Even more striking is the difference regarding the total emission rate of strings of mass m. 
Instead of (|4.17| ) we find 



dFc(m) ^A' D g 2 c Mm- R ^ , (4.21) 



dm 

where = m(M — m)/M and A!' D = 2 ix^ L ^~T n 2 . The power of tur implies that 
light strings are produced much more copiously than heavy ones. However, the factor of M 



12 In order to have the same Haggedorn temperature for open and closed strings we keep a' = 1/2. 
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guarantees that the lifetime is proportional to M 1 , as for open strings. Similarly, the mean 
kinetic energy released per decay is still given by ( [4.151 ). 



As in the case of open strings, there is a simple semiclassical interpretation for ( 4.21| ), 
at least in the limit 1 « m « M. In order for the initial string to split, one point 1 has 
to touch a second point 2. The point 1 can be anywhere along the length L of the string, 
and this explains the factor of M ~ L. Between 1 and 2 there is a piece of string of length 



I ~ m, and if we assume that the string shape is that of a random walk p2|,|18| , the mean 
distance between the two points will be ~ m 1 / 2 . The probability for the two points to meet 
is inversely proportional the corresponding volume 

m (D-l)/2 in 

D — 1 space dimensions. For 



m <C M, mp> « m, and this explains the mass dependence in ( [4.2 



C. Decays of light strings 

By "light strings" we actually mean those with M 2 = 2(N — 1) not large enough to 
use asymptotic approximations for the decay rates.0 Then some of the results obtained for 
heavy strings may no longer be valid. Such is the case with the energy distributions ( [4.9| ). 
( |4.13| ) and (|4.20|) , and with the production rates ( |4.17|) and ( |4.21|) . However, other formulas 



do not rely on the asymptotic approximation ( [4.5|) for the mass degeneracies and may still 
be valid for N in the tens or hundreds, depending on the level of precission sought. This is 
true of eqs. O and ( pLlOD . 

On the other hand, for small values of N one can obtain exact results by summing the 
whole series in ( [4.1|) . A simple example which serves to illustrate the fact, mentioned in 
Section 2, that F(N, N') counts "longitudinal" states for D < 26 is the decay 2 — > 1 + 1 
of the symmetric rank-two tensor ( |2.29| ) into two photons. For photon production by light 



strings, one simply uses ( |4.1| ) with F^u = 0, getting 



13 For simplcity we will consider only light open strings. The generalization to closed strings is 
trivial. 
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F(2, 1) = 0(1) + £(0) = (D - 2) + 1 = D - 1 . 



(4.22) 



This case is so simple that we can compute the amplitudes by hand and sum the probabilities 
over all the polarizations of the tensor and remaining photon according to the definition fl2.1p 
of F(N, N'). We find 

26 - D 



D — 1 — 



8(D-1) 



(4.23) 



which agrees with ( 4.22 ) only for D = 26. The reason for the discrepancy is the following. At 
the critical dimension the only propagating physical state at mass level two is the rank-two 
symmmetric tensor. However, the "longitudinal" state [H| 



r D-l D+l, 

[a-i ■ a_i H — v ■ 0-2 H ^ ' a ~^' 



-p 2 = m 2 = 2 



(4.24) 



satisfies the Virasoro conditions and has positive norm for D < 26 

2 



(010) 
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(D-l)(26-D) . 



(4.25) 



The decay of this state into two photons gives a contribution (26 — D)/8(D — 1) which exactly 
cancels the discrepancy. This example shows that, for D < 26, the double sum F(N, N') 
includes states which are not simply obtained by truncation of the D = 26 spectrum. 

For emission of states with m 2 > we must determine T v {w) by direct expansion of the 
two point-function in powers of v and w. Fortunatelly, we only need to compute a few terms 
in the expansion. The reason is that the degeneracies vanish for negative arguments, and 
by O and ( gl^ the non-universal contribution takes the form 

1 



F NU {N, N')=a u Q N- v{A + 1) - ^m 2 {A 2 + A) 



+ ... 



(4.26) 



One can easily show that 



max 



1 i vn^ 1 

N — v{A + 1) — -m 2 {A 2 + A) ]= 1 + Int(— ) 
2^8 



(4.27) 



so that T v {w) must be known only throught order O ( W u + 1 + Int ( m2 ^/ 8 )) . One can also show that 
the limit ( |4.27|) is approached only for some "soft" decays with very little kinetic energy, i.e. 
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with M ~ m + M', and that for most decay modes the non-universal contribution vanishes 
because the lowest power of w in T u (w) w uA +2 m2 ( A2 + A ) i s greater than N. 

Since the leading universal term in ( |4.1| ) is proportional to Q(N') = Q(l + M /2 /2), 
whereas by ( [4.27]) the leading non-universal term is, at best, proportional to Q(l + m 2 /8), 



non- universal terms are important only for soft decays with m ~ 2M' ~ 2M/3 or greater. 
Indeed, for m > 2M/3, A = in flO|) an d ^(-W, W) = F NU (N, TV'). This is not surprising, 
since for m > 2M/3 the "emitted string" is more appropriately described as a "remnant", 
and the rate is expected to depend on its particular state. 

Most decay modes are accurately described by the universal series in ( |4.1|) and, since 
this is a rapidly decreasing series, it is often enough to keep the first few terms. Moreover, 
there are even decay modes for which only the first term in the series is non-zero. For 
M 2 > 2 (m 2 + M' 2 + 1) the argument of Q for the second term in ( f4.1|) is already negative 
and the following formula is exact 

F(N, TV') = (N-N' - -m 2 ) Q(N') . (4.28) 

As an illustration of these general remarks we will consider the production of m 2 = 2 
states, i.e. decays of the form iV — > 2 + N'. Here we will work at the critical dimension 
D = 26. In this case there are non-universal contributions and we must first compute Xq(w) 
and I\(w) from ( [2.30|) . The result is 



l (w) = (1 + 8w)(^ P n 2 - IQw foX W) + 4w (^D + 0( 
X 1 {w) = -2(w + Qw 2 ) {t^tff + <w + Qw 2 ) (^ P C W) - 4^ 2 (^^) + 0(w 3 ) (4.29) 



where, following the analysis around eq. (|4.27|) , we have included terms of order w and 



w 2 



respectively. The values of the Lorentz invariants for different polarizations of the emitted 
tensor are given in Table I, where the components are referred to the emitted tensor rest 
frame, 1 is the direction in which the tensor is emitted, and a, b = 2, 3, . . . , 25. 
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Polarization 








Cab 








1 

2 


Cla 





1- 


1 

2 




12 T 2 
25 


12 „ 
25 X 


1 

2 



TABLE I. Values to be used in ( ggg ) with x = (1 + n) 2 /2 - M 2 . 



For example, for 7 — > 2 + 3 we have n = 7 — 3 = 4 and ( |3.14j ) gives 



X 4 (w) = 3w 4 + w b + 



which, using Table I with x = ~ yields 



F(7,3) =3^(3) + ^(l) + ^ 



2 = 9626 for £ ab 
= 9624 for f ol 
2 = 9626 for • 

The other possible non-tachyonic decay modes of the form 7 — > 2 + N' are 7 
7 — > 2 + 1. The corresponding sums are universal and given by 

J 5 (w) = 4w 5 + 0(w 8 ) , F(7,2) = AG (2) = 1296 
J 6 (w) = 5w 6 + 0(w 10 ) , F(7,l) = 5G(1) = 120 . 



(4.30) 



(4.31) 



2 + 2 and 



(4.32) 



Note that out of the three decay modes only 7^2 + 3, which is the one with the smallest 
kinetic energy, has non-universal contributions and that these are very small, of the order of 
0.02%. A similar computation for 3^2 + 1 and 4 — > 2 + 1 yields non-universal contributions 
of 8% and 4% respectively. This confirms that these contributions are significant only when 
the emitted particle carries a large fraction of the original mass. Also note that the relative 
contribution of the second term in ( f4.31|) is rather small, of the order of 0.25%. Similar 



examples show that (|4.10| ), which uses only the first term in the solution, gives very accurate 
results for iV as low as 10 or 20. 
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V. DISCUSSION 



In this paper we have considered the problem of computing the averaged semi-inclusive 
rates for emission of specific string states of mass m. This requires the evaluation of the 
double sum 

F(N, N') = J2 E \(®f\V(k)\$i)\ 2 , (5.1) 

where V(k) is the vertex operator for a string state of mass m. Our solution involves a 
recursion relation ( |3.9| ) which is a linear difference equation and depends only on the mass 
m of the emitted string, with the details of the specific emitted state entering only as "initial 
conditions" . The universal series in the solution ( |3.14| ) 

A 

(n — m 2 (r> — 



Z*W = E(n- m\p - 1 )) w ^-W(p 2 -p) + . . . (5.2) 



P =i 

suggests some sort of "deformation" of Plank's distribution for massless strings 

InH = ^- = En^. (5.3) 
1 w p=i 

It would be interesting to check if this series is truly universal and appears in other aver- 
aged semi-inclusive quantities such as iV-point functions, or is just particular to the 3-point 
functions considered here. 

Being able to compute averaged semi-inclusive 4-point functions would also be interesting 
in its own right. They would give us the "form factors" of typical fundamental strings, with 
direct information on their shape and size. This would be an alternative approach, along 
the lines of [23], to the one pursued in [18] for the determination of the sizes of fundamental 



strings [24|. One could also try to compare the form factors with the scattering properties 



of black holes, and see it some of them can be adscribed to a "decoherence" process Ul6| . 

Our analysis has shown that the decoherence procedure of averarging over all initial 
states of mass M is very efficient in suppressing the dependence on the state of the emitted 
string, unless this string carries a very substantial fraction of the original mass. In fact, 
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most decay modes are well described by the first term in the universal series. Dividing the 
approximate expression ( (4.8| ) for F(N,N') by G(N)Q(N') yields the "averaged interaction 
rate" 

< \(<S> f \V(k)\<S> t }\ 2 >^ (N-N'- ^m 2 )-^- = |^ (5.4) 

for the reaction M — > m + M' with p = k + p'. Note that in this case we have averaged over 
the t/iree strings involved. For closed strings the result is similar 

<\(^f\V(k)m\ 2 >^ { -^0- (5.5) 

and was obtained in [|K|] by a different method. These results are very simple, and suggest 
that similar simplifications might take place for higher order interactions. 

Our results for the production rates of strings of mass m by heavy open and closed strings 

dT (m) .. 9 dTJm) ... -£=i m(M-m) , . 

have very simple semiclassical interpretations and provide a microscopic basis for the type 
of heuristic interactions used in the Boltzmann equation approach to hot string gases HQ . 
In particular, the small value of the mean kinetic energy obtained in ( |4 . 1 5|) gives support 
to the assumption that the total length of the strings is approximately conserved. It is 
interesting that the mass dependences in (|5.6| ), which from a semiclassical viewpoint are a 



consequence of the random walk structure of long strings, arise microscopically as a result 
of the interplay between mass degeneracies and phase space factors. Their effects cancel 
out for open strings, giving rise to a flat spectrum. The decay rates for typical open strings 
computed in this paper agree, within a factor of 2, with those computed for strings in very 
special states |||J. 

For open strings on a Dp-brane the spectrum is no longer flat. The string can oscillate 
in D dimensions, but its endpoints are bound to a p + 1 manifold. The phase space factor 
in (4.1U) becomes k p ~ 2 } and instead of a flat spectrum we find a rate proportional to m^ 2 , 



where d = D — p — lis the number of "extra dimensions" orthogonal to the brane. This 
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may be relevant to the decays of "string balls" in brane- world scenarios |19| . It would be 
interesting to study the decay properties of strings in other situations involving D-branes, 
such as those considered in mj. Another possible extensions along these lines would be 
to include the effect of external fields and non-conmutative limits of the theory [[7] on the 
decay rates. 



ACKNOWLEDGMENTS 

It is a pleasure to thank J.M. Aguirregabiria, I.L. Egusquiza, R. Emparan, Jaume Gomis, 
Joaquim Gomis, M.A. Valle-Basagoiti and M.A. Vazquez- Mozo for useful and interesting 
discussions. This work has been supported in part by the Spanish Science Ministry under 
Grant AEN99-0315 and by a University of the Basque Country Grant UPV-EHU-063.310- 
EB187/98. 

APPENDIX A: 

Here we will compute the integral ( ^.8[ ) which appears in the derivation of the recursion 
relation. The scalar two-point function on the cylinder ip(v,w) is related to the Jacobi $1 
function by [13| 



4(v,w) = — -±±J- , T= — , , v=- (Al) 

t (U|7~j mw mw 

and using the product formula for the theta function yields the following expression 

2vr -^-l-2g 2n cos27rz/ + g 4n ,, . 

^,w) = -—mncu]l (1 _ q2n)2 , (A2) 

where q = e lnT . The function tp(v, w) has zeroes for v = w n with n — 0, ±1, ±2, . . ., and 
expanding (|A2]) around v = 1 gives 

tjj( V} w ) = - \n(v) + 0((ln v) 3 ) = -r + 0(r 3 ) , (A3) 
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where r is the euclidean proper time used throughout the paper and should not be confused 
with the modular parameter f in (|AI|). This fixes the behavior of the function ip(v, w) near 
v = 1 

f-\n 2 v\ 1 

ifi(v, w) = ^j exp — ifi(v, w) = -r(l + -r) + 0(r 3 ) (A4) 

\2mw J 2 

and implies 



-*-Mi) e *-*«.W) = ^ w )™ 2 = r - 2 (i + + 0(r 2 )) 
(<9 T ^(l)c> T X») = rr%fo$(v,w) = ~Q(v,w) = -1 + 0(1) . (A5) 



e 



T 2 



Each additional proper time derivative contributes a negative power of r, and (|A5|) together 
with the mass shell condition implies the following form for the correlator 



A 2 

(V'\k, 1) V'(k, v)) = -(1 + ^-r) + 0(1) . (A6) 



The value of A can be found by noting that computing the leading singularity of the 
correlator with the help of ( |A5| ) is equivalent to finding the leading term in the OPE of 



ordinary (unprimed) vertex operators. But this term is given by 

K t m-(r) = ^ + ... (A7) 

where is essentially Zamolodchikov's inner product, which is related to the familiar 

quantum mechanical inner product^ by = Comparing (|A6|) and (|A7]) implies 

A = 1 for normalized states. Computing the integral gives 

r du j- r im 2 m 2 

A — V "(y'\k, 1) V'{k, v)) = f dre nT -(l + —r) = n + — . (A8) 

JC V JO T A 1 1 
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See chapter 6 of [14| for details 
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APPENDIX B: 



In this Appendix we present an approximate expression for the mass level degeneracy 
which is more accurate than the usual asymptotic formula (|4.5|) . We have used it in our 
numerical computations in Section 5, and find it particularly useful for the range of values 
of N where a direct power expansion of the partition function is not practical, while the 
asymptotic formula is not yet accurate enough. 

The basic idea is to compute the first corrections to the saddle-point aproximation for 

r din 00 
G(N, D) = i -^ W - N f{ w f- D , f( w ) = n (1 - «,») . (Bl) 

J W n=l 

Since the saddle-point is close to w — 1, one has to use the modular transformation 

/H=( ] ^) 1/ V 1 / 2 V /12 /(g 2 ) , g = e 2 ^ (B2) 

which gives 

ln/H = ~ hi2vr - \ In/? + 1/J - ^ + 0{e~^^) (B3) 

where we have made the change of variable w = . Keeping all corrections of order 
0(N~ l l 2 ) yields the following result 



where 



r x / £>-2 x l D-2 f , , \ 

a(x) = (N + ( 4vr 2 x - 12 In x - 12 In 2tt ) 

K J K 24 x 24 V y 

6(x) =4x 2 (2tt 2 x-3) 

9 4vr 2 x- 3 45(tt 2 x-1) 2 
C[X) = 2(2tt 2 x-3) 2 ~ (2tt 2 x-3) 3 ( } 



and the position of the saddle-point is given by 
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The following table gives a few sample values for the degeneracies, and we can compare 
the exact results for Q(N,D) with those obtained from (|B4j) . Note that the agreement is 
very good for N as low as 20, and even the results for N — 1 are acceptable. 



N 


G(N, 26) 




G(N, 16) 


( 


B4 


) 


1 

5 

20 
100 


24 
176,256 
2.16109- 10 14 
2.43478 • 10 40 


23.819 
176,271 
2.16118- 10 14 
2.43480 ■ 10 40 


14 
19, 754 
3.41780- 10 11 
7.90265 • 10 31 


13.949 
19, 759 
3.41806 • 10 11 
7.90276 ■ 10 31 



TABLE II. Sample values of G(N, D) 



This should be compared with the usual asymptotic approximation 

G(M) ~ (27rT H )- ^M-^e M ' T » , T H = -J^- (B7) 
which, for = 100, gives 4.79 • 10 40 for D = 26 and 9.20 • 10 31 for D = 16. 
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